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Coherent radiation at radio frequencies from high-energy showers fully contained in a dense radio¬ 
transparent medium - like ice, salt, soil or regolith - has been extensively investigated as a promising 
technique to search for ultra-high energy (UHE) neutrinos. Additional emission in the form of tran¬ 
sition radiation may occur when a neutrino-induced shower produced close to the Earth surface 
emerges from the ground into atmospheric air. We present the first detailed evaluation of transition 
radiation from high-energy showers crossing the boundary between two different media. We found 
that transition radiation is sizable over a wide solid angle and coherent up to ~ 1 GHz. These prop¬ 
erties encourage further work to evaluate the potential of a large-aperture UHE neutrino experiment 
based on detection of transition radiation. 


I. INTRODUCTION 

The nature and origin of UHE cosmic rays (UHECRs) 
with energies around 10^° eV is one of the most puzzling 
questions in particle astrophysics mm- Charged cosmic 
rays with energy less than ^ 10 EeV are significantly 
scattered by the galactic magnetic held and cannot be 
traced back to their origin. At higher energies where de¬ 
flections become less important, the cosmic ray flux falls 
steeply Em, requiring giant experiments like the Pierre 
Auger Observatory [5] and the Telescope Array [6] to 
collect reasonable statistics. Neutrinos, which travel un- 
dehected through the universe, may provide important 
clues on the origin of UHECRs [7]. In fact, cosmic rays in¬ 
teracting with matter and radiation present at the accel¬ 
erating source, or encountered along their path towards 
Earth, may produce UHE neutrinos. Recently, extrater¬ 
restrial neutrinos of energies up to few PeV have been 
detected by the IceCube Neutrino Observatory [8j. For 
neutrino energies in the EeV range, a promising detection 
technique is based on the radio signal from the neutrino- 
induced shower in a dense radio-transparent medium. Ice 
has been proven to be an ideal medium for this tech¬ 
nique by pioneering experiments such as RICE [5] , which 
buried radio antennas in the South Pole ice, and ANITA 
[lOj . which looks with balloon-flown radio antennas over 
a large portion of the Antarctica ice cap. Large ground 
arrays of antennas, buried in ice (ARA [TT]) or on the ice 
surface (ARIANNA [12] and GNO [12]) are now under 
development. In another interesting approach, several 
groups [UHE] have pointed radiotelescopes towards the 
Moon, searching for radiopulses from neutrino interac¬ 
tions in the lunar regolith. All of these experiments are 
based on the Askaryan effect m - the emission of co¬ 
herent Cherenkov radiation from the ~ 20% excess of 
electrons present in an electromagnetic shower. 

In this paper we explore transition radiation (TR) [IS] 
at radio frequencies as a possible way to detect UHE neu¬ 
trinos. Transition radiation is emitted by charged parti¬ 
cles when crossing the boundary between two media with 


different indices of refraction. From the theoretical view¬ 
point, TR can be understood as the radiation required 
to match at the boundary the electric field produced by 
the charged particle in the bulk of the two media. It has 
been well studied both theoretically and experimentally; 
for a comprehensive review of transition radiation, see for 
example |19j . However, there are only a few estimates in 
the literature of TR from cosmic ray showers - either 
crossing from atmospheric air to a denser medium (e.g. 
clouds |20j or ground |2Ij l or emerging from ground into 
air (a rough estimate is given in [12]). While the meth¬ 
ods presented in this paper are of general application, 
we are mostly motivated by the study of TR emitted 
when a shower, originated by a UHE neutrino interac¬ 
tion below the Earth surface, escapes the ground into 
the atmosphere. In fact, the net electric charge corre¬ 
sponding to the excess electrons of the shower will pro¬ 
duce an upward-going transition radiation when passing 
through the ground - atmosphere boundary. Since show¬ 
ers in dense media are rather compact in size, we expect 
the TR to be coherent up to GHz frequencies, with signif¬ 
icant advantages in the cost of detectors and background 
noise level. In this work we present the first detailed eval¬ 
uation of transition radiation from high-energy showers 
crossing the boundary between two different media. 

In Section [TT] we present a general method to calculate 
transition radiation from a shower crossing the boundary 
between two media using the far-held limit. This method 
extends the well-known Zas-Halzen-Stanev (ZHS) algo¬ 
rithm [23] to account for boundary matching of the ra¬ 
diation held. In ZHS, the trajectory of each particle in 
a homogeneous medium is approximated by linear seg¬ 
ments, sufficiently small that the particle can be assumed 
to move uniformly, and the corresponding electromag¬ 
netic radiation is computed exactly in the Fraunhofer 
limit from the Maxwell equations [23]. In the extended 
ZHS algorithm, hereafter called ZHS-TR, a particle seg¬ 
ment crossing the two media is split at the boundary [22] 
and rehection and refraction of the radiation emitted by 
the two track segments is then included in a standard 
way |2S]. 
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In Section III the properties of TR emitted by a high 
energy shower when crossing from a dense medium to 
air are studied in detail. In particular, we investigate 
the spectral characteristics of the TR signal and their 
dependence on the shower energy and zenith angle, on 
the stage of the shower development at the boundary, on 
the energy of the shower particles, and on the type of 
dense medium. 

Conclusions are drawn in Section IIVI 


II. TRANSITION RADIATION FROM 
HIGH-ENERGY SHOWERS: CALCULATIONS 


superposition of the electric field from each sub-track, 
with relative phase shifts due to different starting point 
positions and time delays properly taken into account. 

The electric field E{t,x) produced by a charged par¬ 
ticle moving with uniform velocity v between two fixed 
points in a homogeneous medium can be derived from 
Maxwell equations, and has a frequency Fourier trans¬ 
form [32] 


E{ui,x) = 2 / E(t,x)e^‘^*dt 


( 1 ) 


given by [23] : 


A. The ZHS Monte Carlo 

Our calculation of radio emission from particle show¬ 
ers crossing the boundary between two media is based on 
the methods ^31 [331 [37] implemented in the ZHS Monte 
Carlo, which simulates electromagnetic showers and their 
associated coherent radio emission up to EeV energies 
[28j . Originally developed for the Fraunhofer limit in 
homogeneous ice [35], it has been extended to other di¬ 
electric homogeneous media [35] [SO] and to reproduce 
near field effects by dividing the particle trajectories in 
sufficiently small sub-tracks [3T]. To provide the reader 
with the necessary background, the methods currently 
implemented in the ZHS Monte Carlo are reviewed in 
the following. 

The ZHS Monte Carlo follows electrons and positrons 
interactions down to a kinetic energy threshold of 100 
keV. Charged particles below this threshold contribute 
negligibly to the electric field, which is proportional to 
the particle tracklength [33]. The simulation includes 
bremsstrahlung, pair production, and the interactions re¬ 
sponsible for the generation of the excess charge (Mpller, 
Bhabha, Compton scattering and electron-positron an¬ 
nihilation). In addition, multiple elastic scattering (ac¬ 
cording to Moliere’s theory) and continuous ionization 
losses are implemented. The track segment between two 
consecutive particle interactions is divided into multiple 
linear segments, along which the particle is assumed to 
move with constant velocity. To ensure a proper calcu¬ 
lation in the Fraunhofer limit, the maximum equivalent 
depth of each sub-track does not exceed 0.1 radiation 
lengths. Also, for low energy particles the size of the sub¬ 
track is required to be smaller than the particle range, 
and this step size is used to evaluate ionization losses and 
multiple elastic scattering. Accurate bookkeeping of the 
absolute timing is maintained during particle propaga¬ 
tion, including geometrical time delays as well as those 
introduced by different particle velocities (a uniform en¬ 
ergy loss along the sub-track is assumed for this purpose). 
An approximate account is also made of the time delay 
induced by multiple elastic scattering. From the position 
and time of the endpoints of the sub-tracks, the frequency 
spectrum of the electric field is derived. The total elec¬ 
tric field at the observer’s location is then calculated by 


E{u;,x) = 


enr 


27reoc^ 


pikR 

R 




i{(jj — k ■ v) 

( 2 ) 


In Eq. ([^, i? is the distance from the track to the 
observer, which is assumed to be large enough for the 
Fraunhofer regime to be valid. Wave vector k points in 
the direction from the track to the observer and has mod¬ 
ulus k = Lu I o', with c' the speed of light in the medium. 
Also, v± = —kx (kxv) is the component of v in the plane 
perpendicular to the unit vector k — k/\k\. Finally, ti 
and t 2 = ti + St are the absolute times at which the par¬ 
ticle passes through the starting point and endpoint of 
the sub-track, respectively. In magnetic materials, fir is 
the relative permeability of the medium. 


B. The ZHS-TR algorithm 

1. General case: Transition between two media 


The original ZHS Monte Carlo described in Sect. HA 
has been mostly used in simulations where the shower is 
contained in a single dense medium, which also hosts the 
observer. Several modifications are required to simulate 
a shower extending over two different media. In the fol¬ 
lowing, we will detail the ZHS-TR algorithm for a shower 
starting in a medium of refractive index rii, separated by 
a planar boundary from a second medium of refractive 
index n 2 where the observer is located. This configu¬ 
ration accounts for the case of a shower initiated below 
the Earth surface by an UHE neutrino interaction, and 
emerging in the atmosphere where it is detected. How¬ 
ever, the algorithm is general and can be used for other 
cases (for example, the shower could start in air, or the 
observer could be located in medium 1). 

As in the original ZHS algorithm, a particle track is ap¬ 
proximated by linear segments along which the particle 
is assumed to move with constant velocity. However, the 
electric field from a given sub-track can include one, two 
or three contributions, depending on whether the sub¬ 
track is fully contained in medium 1, in medium 2 (where 
the observer is) or crossing the boundary. The different 
contributions are pictorially represented in Fig. [^ We 
calculate the electric field at the observer position, with 
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FIG. 1. Schematic illustration of the ZHS-TR algorithm (see 
text for a detailed description). The three cases shown corre¬ 
spond to a sub-track crossing the interface, a snb-track in the 
same medium as the observer (track above) and a sub-track 
in the other medium (track below). 


the sub-track viewed in the direction k. Also, we assume 
that Fraunhofer conditions apply, namely XR ^ and 
R ^ D, where A is the wavelength of the radiation and 
D the characteristic size of the sub-track. With this ap¬ 
proximation the factor 1/i? in Eq. ([^ can be factored 
out as if it was constant. 

For tracks fully contained in medium 1, the only con¬ 
tribution comes from the electric field (Eq. (§) emitted 

into direction k" (the direction of a ray refracted at the 
boundary into direction k according to Snell’s law, see 
Eig. [^. The radiated electric field is then decomposed 
into s and p polarizations (perpendicular and parallel 
to the plane of incidence, respectively), and transmitted 
through the boundary applying the appropriate Fresnel 
coefficients [33]. Because we are dealing with a point 
source and not plane waves, the standard transmission 
Fresnel coefficients must be multiplied by a factor 


n 2 cos 9 


Til / 


2 

r- 

^sin6» 

ni 



(3) 


which can be understood as describing the change in di¬ 
vergence of rays upon transmission |24j . The angle of 
refraction 6 is the angle between the unit vector k and 
the normal to the boundary plane. The final transmis¬ 
sion formulas read 
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When calculating the electric field at the observer posi¬ 
tion, the p polarization in the reflection plane remains of 


course perpendicular to k, the new propagation direction 
after refraction. In the following we refer to the electric 
field corresponding to tracks fully contained in medium 1 
as the “transmitted contribution”. 

For tracks fully contained in medium 2, two distinct 
contributions must be evaluated. The first is the stan¬ 
dard contribution given by Eq. describing the elec¬ 
tric field radiated into direction k (“direct contribution”). 
The second contribution (“reflected contribution”) comes 
from radiation emitted into a direction k which is re¬ 
flected off the boundary into the observing direction k 
(Fig.0. Eq. ([^ is also used to calculate this contribu¬ 
tion. The electric field is again decomposed into s and 
p polarizations, and the appropriate Eresnel coefficients 
are used for reflection: 


ni cos 0 — 712 
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(5) 


Since reflection does not change the divergence of rays, 
Eqs. ([^ have their usual plane-wave form. As in the 
transmitted contribution, the p polarization remains per¬ 
pendicular to k after reflection. 

Finally, for tracks crossing the two media the following 
procedure is applied. The track is split at the boundary 
plane into two sub-tracks, as in [33]. Each sub-track is 
then fully contained in one medium and the methods pre¬ 
viously described can be applied to calculate the electric 
field. For these tracks transition radiation is naturally 
accounted for with this procedure. In the ZHS algorithm 
each of the endpoints of a sub-track contributes a term 
to the electric field. For a charged particle moving uni¬ 
formly in a homogeneous medium (i.e. with constant 
velocity along the track), the contributions from the com¬ 
mon endpoints of two adjacent sub-tracks cancel exactly. 
Thus, the electric field associated to the track does not 
depend on the number or on the length of its sub-tracks. 
On the other hand, a track crossing two media has two 
adjacent sub-tracks with common endpoints located at 
the boundary. Since these endpoints are associated to 
sub-tracks contained in media with different index of re¬ 
fraction, their contributions do not cancel exactly. Tran¬ 
sition radiation appears because of this incomplete can¬ 
cellation. Notice that the assumption of uniform velocity 
can be made to an arbitrary degree of precision by reduc¬ 
ing the length of the sub-tracks. 

The electric field produced by the entire shower is then 
obtained from the superposition of the individual contri¬ 
butions of all particle tracks. Explicit geometric phase 
differences between tracks on both sides of the bound¬ 
ary plane are appropriately taken into account in this 
procedure. This concludes the ZHS-TR algorithm. 
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The technical implementation of the ZHS-TR algo¬ 
rithm requires two simulation runs, since the origi¬ 
nal ZHS code was developed for a single homogeneous 
medium. The simulation starts by generating and prop¬ 
agating the shower in medium 1. Then, particles crossing 
the boundary are propagated in medium 2 in a separate 
simulation. The ZHS code presents some limitations, 
since it does not treat hadronic interactions and radio 
emission induced by the Earth magnetic held |H4j . and 
assumes a constant density of the media, which is not ap¬ 
propriate for atmospheric air. However, neglecting these 
effects should not impact signihcantly our results as dis¬ 
cussed in Sect. nni 



2. Vacuum approximation 


In the case under consideration, transition radiation 
is produced when the shower exits into the atmosphere. 
Since the air density is much lower than the density of 
the medium where the shower originated, medium 2 can 
be reasonably approximated by vacuum to simplify cal¬ 
culations. There is no difference in the simulation for 
the portion of the shower contained in medium 1: parti¬ 


cle tracks are treated as in Sect. HB 1 and their corre¬ 


sponding transmitted contribution is calculated. Instead, 
particles crossing the boundary, which do not interact in 
vacuum, are modeled as moving to infinity with constant 
velocity. The electric field associated to such a semi¬ 
infinite track [24] in the Fraunhofer regime is given by: 

pikR ,, i(ui-k-v)ti 

1-1/ \ C-fX'f . C U J_C. , , 

E{uj,x) = - -luj - 5 ---—, ( 6 ) 

27 reoc 2 R i{uj — k ■ v) 

where the notation follows Eq. ([^. The origin of this 
formula can be understood by dividing the semi-infinite 
track into sub-tracks. The electric held is obtained by 
summing Eq. ([^ over all sub-tracks. Since contributions 
of the common endpoints of adjacent sub-tracks cancel, 
the only term left corresponds to the starting point (the 
term at inhnity goes as 1/i? —0), thus giving Eq. (§. 

Notice that since the electric held of the entire track is 
given by Eq. ([^, there is no need for a division in sub¬ 
tracks, which makes the simulation faster. Equation (§ 
is then used to calculate the direct and rehected contri¬ 
butions to the electric held following the procedure of 
Sect. 


Sect. 


HBl 

Tim 


We will apply the vacuum approximation in 


C. Transition radiation from a single particle 


FIG. 2. Angular distribution of the far-held electric held 
radiated by a single charged particle crossing from ice to air: 
comparison of the results obtained using the ZHS-TR algo¬ 
rithm with sub-tracks 10 cm (crosses) and 100 cm (open dots) 
long to the exact calculation (black line). 


0.99c, corresponding to an energy of about 3.6 MeV. In 
order to approximate an inhnite track (required for the 
classical derivation), electron interactions were not con¬ 
sidered, and a total track length of 100 km was equally 
divided between the two media. We then evaluated the 
electric held at an observer position d of zenith angle 0, 
and |d| = 1 km (the origin of the coordinate system coin¬ 
cides with the intersection of the track with the boundary 
plane). Since this distance is much larger than the length 
of the sub-tracks in ZHS-TR, we expect the Fraunhofer 
limit used in Eq. © to be a good approximation. 

In Fig. 1^ the electric held calculated with our algo¬ 
rithm is compared with the exact result obtained by 
solving the Maxwell equations for an inhnite track |25j . 
The two approaches are in practically perfect agreement 
when sub-tracks of 10 cm length are used in ZHS-TR. 
To illustrate the importance of the Fraunhofer limit, we 
also show the same comparison for sub-tracks of 100 cm 
length. At small zenith angles, where a discrepancy is 
observed, the perpendicular distance between the track 
and the observer (~ |d|0) becomes comparable to the 
length of the sub-tracks, and the Fraunhofer regime is no 
longer valid. 

In the rest of the paper, we will assume that the ob¬ 
server is always distant enough from the shower for the 
Fraunhofer regime to be valid. Due to the finite extent 
of the shower, this is a good approximation in most prac¬ 
tical cases, for example in the detection of the TR signal 
from a satellite or balloon instrument [TO], [35]. 


A basic check of the ZHS-TR algorithm was per¬ 
formed by estimating the transition radiation from a sin¬ 
gle charged particle and comparing the result with the 
classical derivation [25]. An electron crossing the bound¬ 
ary from ice to air in the vertical direction was consid¬ 
ered for this purpose. The electron velocity was fixed to 


III. PROPERTIES OF TRANSITION 
RADIATION FROM HIGH-ENERGY SHOWERS 

In this section, we investigate in detail the properties 
of the radiation emitted by an electromagnetic shower 
transitioning from a dense medium to air. Showers ini- 
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tiated by a 100 TeV electron in ice (ni = 1.78 at radio 
frequencies) were used in most of these studies. Also, no 
thinning was employed in the simulation of the showers 
{i.e. all particles were propagated). Since in the Fraun¬ 
hofer approximation the electric field E(a;,x) is inversely 
proportional to the distance R to the shower, results are 
presented in terms of the quantity \R x E|. The case of 
a vertical shower is treated first, and then results for in¬ 
clined showers - which may be more realistic for a UHE 
neutrino-induced interaction - are presented. 


A. Vertical showers 


1. Direct, reflected and transmitted contributions 


The angular distribution of the electric field radiated 
by a 100 TeV vertical shower simulated with ZHS-TR is 
shown in Fig. (black solid line). The shower’s start¬ 
ing point in ice was chosen appropriately for the shower 
maximum to occur approximately at the ice-air bound¬ 
ary. The zenith angle dependence reflects the contribu¬ 
tions to TR by single electrons whose angular features 
are shown in Fig. The weak angular dependence of 
the emitted radiation is a notable feature, when com¬ 
pared to coherent Cherenkov radiation which is strongly 
beamed around the Cherenkov angle |23]. A value of 
|i? X E| ^ 5 X 10“^V/MHz is obtained for this shower. 


Also shown in Fig. are the contributions to the to¬ 
tal electric field (direct, transmitted and reflected) intro¬ 
duced in Sect. Em The direct contribution from parti¬ 
cle tracks in air is dominant at all angles, and presents a 
peak at ~ 2°. This peak can be understood as Cherenkov 
radiation produced in air (Cherenkov angle ~ 1.4°) by 
the high energy shower particles which follow closely the 
shower axis. In addition, there are noticeable fluctua¬ 
tions around the mean value of |i? x E|, which are due 
to the fine structure of the shower development in air. 
In fact, the spatial dimension corresponding to the fre¬ 
quencies of observation - 200 MHz and 1 GHz - is much 
smaller than the transverse size of the shower (a simu¬ 
lation performed at 10 MHz, not shown for the sake of 
brevity, presents a much smoother behavior). The re¬ 
flected contribution is negligible except for the largest 
angles, where it becomes comparable to the direct con¬ 
tribution. However, the two contributions cancel each 
other for an observer looking orthogonally to the shower 
axis, since a tt phase shift is gained upon reflection, and 
the total electric field drops at zenith angles approaching 
90°. Lastly, the transmitted contribution from particles 
in ice is never dominant for a vertical shower. In fact, 
the coherent Cherenkov radiation from the Askaryan ef¬ 
fect in ice undergoes total internal reflection at the ice- 
air boundary plane (the Cherenkov angle is larger that 
the critical angle for ni=1.78) and cannot reach the ob¬ 
server. However, this component can become significant 


for inclined showers, see Sect. HIB 


8. Particles crossing the boundary and vacuum 
approximation 


To further investigate the origin of the emitted ra¬ 
diation, we performed a dedicated ZHS-TR simulation 
where the electric field was calculated taking into ac¬ 
count only contributions from particles crossing the ice - 
air boundary. Results are shown in Fig.|^ where we plot 
the ratio of the electric field magnitudes obtained in this 
simplied simulation and the electric field magnitudes of 
the full result. Indeed, transition radiation from particles 
crossing the boundary accounts for most of the emission 
(at low zenith angles, Cherenkov radiation from particles 
propagating in air is dominant, as seen in Fig. |^. 

Also shown in Fig. |^is a ratio of the electric field mag¬ 
nitudes resulting from a ZHS-TR simulation using the 


vacuum approximation of Sect. HB2 with the electric 


field magnitudes of the full result. The vacuum approx¬ 
imation reproduces the full ZHS-TR simulation reason¬ 
ably well, with a notable exception at low zenith angles. 
This can be explained by the fact that, for an observer 
looking from the vertical direction, highly relativistic par¬ 
ticles with k \\ V have uj — v ■ k ^ 0, resulting in a large 
contribution from Eq. §. High energy particles emerge 
from ice strongly beamed along the vertical axis, and 
since their direction does not change in the vacuum ap¬ 
proximation, a strong peak is expected. On the other 
hand, the full ZHS-TR simulation includes particle in¬ 
teractions in air which randomize the directions of the 
sub-tracks, leading to a suppression of the peak. In any 
case, the agreement between the vacuum approximation 
and the full simulation over most of the zenith angles 
suggests that the details of the shower development in 
air have a marginal effect on the TR emission. Thus, we 
expect that neglecting the altitude dependence of the air 
density and geomagnetic effects in particle propagation 
will not change significantly our estimates of transition 
radiation. 


3. Frequency spectrum 

Another important property of the radiation is its fre¬ 
quency spectrum, shown in Fig. for 100 TeV show¬ 
ers. To reduce shower-to-shower fluctuations, a hundred 
showers were simulated, and the modulus of their sig¬ 
nal was averaged. All averages will be done this way 
unless explicitly stated. The rise at the highest frequen¬ 
cies is caused by contributions from incoherent emission 
and continues well beyond the frequencies displayed. To 
check that the coherent contribution is indeed quickly 
dying away after certain threshold frequency, we fitted a 
power law dependence to the magnitude of the electric 
field above 20 GHz where the coherent contributions are 
negligible. When we subtracted this contribution from 
the total electric field observed in the simulations, the 
remainder showed the expected steep drop off. 

Unlike coherent Cherenkov radiation, which drops off 
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FIG. 3. Angular distribution of the electric field radiated by a single 100 TeV vertical shower crossing from ice to air: scaled 
magnitudes of the total field (black), direct (green short dashes), reflected (red dot-dashes) and transmitted (blue dashes) 
contributions. Due to phase differences the total magnitude can be smaller than those of individual parts. See the text for 
details. The electric field is evaluated at 200 MHz (left panel) and 1 GHz (right panel). 




FIG. 4. Ratio of the electric field magnitudes calculated in the vacuum approximation and the full simulation (solid) and 
the same comparison for simulation where we include only tracks crossing the ice-air boundary (dashed). The vertical shower 
started with 100 TeV energy and crossed from ice to air. The electric field is evaluated at 200 MHz (left panel) and 1 GHz 
(right panel). 


at relatively low frequencies away from the Cherenkov 
angle |23j , transition radiation is coherent up to ~ 1 GHz 
over a wide range of angles. The physical origin of 
this effect may be understood with a simplified model 
where TR is produced only by particles crossing orthog¬ 
onally the boundary plane with velocity c, uniformly dis¬ 
tributed over a disk of radius tm ~ 11 cm (the Moliere 
radius in ice). Under these approximations, each parti¬ 
cle contributes equally in amplitude, and only the rela¬ 
tive phases due to propagation are important. The prob¬ 
lem then becomes analogous to calculating the diffraction 
pattern of a uniformly-illuminated circular aperture |33j . 
Thus, the electric field at given frequency / goes as 


where Ji is a Bessel function of the first kind. This ex¬ 
pression is maximal for / —>■ 0 and reaches the first zero 
when 


0.61c 
tm sin 9 


( 8 ) 


This equation may be taken as an estimate of the TR cut¬ 
off frequency. Indeed, the trend observed in Fig. with 
lower cutoff frequencies corresponding to larger angles of 
observation, is well described qualitatively by Eq. (|^. 


4- Scaling with number of particles crossing the boundary 


|E| 


c 

fsinO 


Ji 


/ 27r/rM 



(7) 


We have shown in Sect. lIIIA2l that most of the emis¬ 
sion from a shower emerging from ice into air is due to 
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FIG. 5. Frequency spectrum of the electric field radiated by 
a 100 TeV vertical shower crossing from ice to air, at 10° 
(black solid), 20° (red solid), 40° (green dot dashes), 60° (blue 
short dashes), and 80° (black long dashes) zenith angle of 
observation (top to bottom). Each curve is an average of a 
hundred showers. The rise at the highest frequencies comes 
from incoherent emission. 


particles crossing the boundary. These particles should 
contribute coherently to the radio signal up to a cutoff 
frequency related to the geometrical size of the shower. 
Thus, we expect the averaged electric field {R x E) to 
scale linearly with the charge excess AQ (the number of 
electrons minus the number of positrons) flowing through 
the boundary. 

To test this we used an averaged value 
{ExR) = ^[ |i? X £;(200MHz,6»)|dn 

= — / |i? X £;(200MHz,6»)|sin6»d6», (9) 

^ Jio° 

where the normalization constants are 


C' = ^ = 

27r 


r70° 


sin 6d9. 


' 10 ° 


( 10 ) 


This quantity measures an average magnitude of the elec¬ 
tric field in the middle angular region which is not af¬ 
fected by the air Cherenkov emission or the decrease of 
the signal at large zenith angles. 

The theoretically predicted linear scaling is confirmed 
in Fig. which shows the dependence of the average 
magnitude of the electric field as a function of the charge 
excess for a hundred simulated 100 TeV showers starting 
at a fixed distance 500g/cm^ from the ice-air boundary. 
One can notice the presence of showers with charge excess 
of only 60% of the maximum observed value. This effect 
is related to the fact that while for most of the showers 
500g/cm^ corresponds approximately to the distance to 
the shower maximum, for others the shower maximum 
occurs at a somewhat different distance. This affects the 
number of crossing particles and relatedly the size of the 
mean electric field. 



AQ crossing the interface 


FIG. 6. Relation between the absolute charge excess crossing 
the boundary and the averaged electric field for a hundred 
100 TeV showers starting 500 g/cm^ from the ice-air bound¬ 
ary. The straight line is a fit of the expected theoretical de¬ 
pendence. 


On top of the systematic increase and decrease of the 
mean value of the electric field with the charge excess 
through the boundary there is also the effect of the fine 
structure of an individual shower. This leads to fluctua¬ 
tions around the mean value, seen for example in Fig. 
From this figure we can estimate the magnitude of these 
fluctuations for a 100 TeV shower; for showers of higher 
energies these fluctuations are relatively smaller due to 
higher total number of particles in these showers. 

In another related study, we simulated 100 TeV show¬ 
ers starting at different ice depths. Since the number of 
particles crossing the boundary depends in this case on 
the stage of the shower development, the signal should 
be largest when the maximum of the shower development 
occurs close to the boundary for a typical shower. This is 
indeed confirmed in Fig. The peak at low observation 
angles (< 5°), which is due to Cherenkov radiation from 
particles in air, becomes less pronounced the farther the 
shower starting point is from the boundary. In fact, at 
later stages of the shower development the shower is less 
compact and its particles have a larger angular spread, 
resulting in a loss of coherence. 


5. Dependence on shower particles energy spectrum 


Since the electric field depends on the particle veloc¬ 
ity (c./. Eq. ([|)), we may expect a dependence of the 
TR on the energy spectrum of the shower particles. To 
study this dependence, we estimated the electric field for 
100 TeV simulated showers taking into account only the 
contributions of sub-tracks from particles with energy be¬ 
low a given threshold. Results are presented in Fig. 
(note the linear scale of the vertical axis). The highest 
energy particles (> 100 MeV) contribute mostly at low 
zenith angles. In fact, these particles are strongly beamed 
in the vertical direction, and account for the coherent 
Cherenkov peak as discussed in Sect. Ill A 1 The con- 










Observer zenith angle (deg) 


FIG. 7. Angular distribution of the electric field radiated by a 100 TeV vertical shower crossing from ice to air, with starting 
point in ice at 300 (blue long dashes), 400 (blue short dashes), 500 (black, showers typically crossing at shower maximum), 
600 (green dot-dashes), 700 (red short dashes) and 800 (red long dashes) g/cm^ from the ice-air boundary. The electric field 
is evaluated at 200 MHz (left panel) and 1 GHz (right panel). Each curve is an average of a hundred showers. 


tributions from particles in the energy ranges 3-10 MeV, 
10-30 MeV, and 30-100 MeV are comparable. This is 
due to a compensation between the number of particles 
- which is higher at low energies - and the amount of 
radiation coming from a single particle - which is higher 
at high energies. Particles with energies below 3 MeV do 
not contribute significantly. 


6. Dependence on the shower energy 

We expect the absolute charge excess for showers which 
cross the boundary at their maximum developments to 
be a good proxy of the shower energy as the number of 
particles at the shower maximum is roughly proportional 
to the shower energy. As a consequence, we expect the 
electric field to scale linearly with the shower energy. 

In our implementation of the ZHS-TR algorithm it is 
not possible to know a priori whether a given simulated 
shower will hit the boundary around the shower maxi¬ 
mum or whether it will be one of the downward fluctu¬ 
ations seen in the Fig. To quantitatively asses how 
the radiated electric field depends on the shower energy, 
we then compare maximal signal seen in all our simu¬ 
lated showers at three different energies: 1, 10 and 100 
TeV. As a first step we run several trial showers to find 
for each of these energies an approximate starting point 
where we can expect showers reaching the boundary at 
roughly their shower maximum. We then simulate a hun¬ 
dred unthinned showers starting at this point and search 
for the maximal value of {E x R). We take this quan¬ 
tity (denoted as {E x as our proxy for the value of 

{E X R) for a shower hitting the boundary at the shower 
maximum. Similar other choices are possible but lead all 
to comparable results. 

The maximal values of the mean electric field at these 


three energies observed in our sample of showers are 

{E X i?)“oo’TeV = 6-4 X 10"^ V/MHz (11) 
{E X i?)foTeV = 7.8 X 10"® V/MHz 
{E X i?)“TeV = 1-2 X 10"® V/MHz. 

They deviate slightly from the theoretical prediction, 
which we attribute to higher effect of fluctuations at lower 
energies. This viewpoint is supported by Fig. which 
shows the angular distributions of the electric fields av¬ 
eraged over the investigated showers. While the 100 TeV 
showers show profiles comparable to those of a single par¬ 
ticle, at lower energies the fluctuations lead to a change 
in the shape of the profile. Notice that even for 100 TeV 
showers the agreement with the single particle profile is 
worse at higher frequencies due to the lower level of co¬ 
herence. 

As a further study, we investigated several thinned 
1 EeV showers. After finding an approximate shower de¬ 
velopment D to the shower maximum we then repeated 
the simulation for showers with starting point separated 
by distance D from the interface. Due to increased com¬ 
putational demands we simulated only twelve showers for 
this energy. The largest value of {E x R) we found among 
these showers is 

{E X = 1.52mV/MHz « 2400 {E x 

( 12 ) 

The reason why the electric field at very high energies 
does not scale linearly with the energy can be traced 
back to the Landau-Pomeranchuk-Migdal (LPM) effect 
[31 [37]. As was found in [33], due to elongation of the 
shower at EeV and higher energies the number of parti¬ 
cles around shower maximum for a shower with the LPM 
effect included decreases with respect to the same shower 
with no LPM effect. Consequently, the charge excess and 
the radiated electric field are also lowered. 
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FIG. 8. Angular distribution of the electric field radiated by a 100 TeV vertical shower crossing from ice to air (from top to 
bottom): all particles (black), particles with energy <100 MeV (blue short dashes), 30 MeV (red), 10 MeV (green dot dashes) 
and 3 MeV (black long dashes). The electric field is evaluated at 200 MHz (left panel) and 1 GHz (right panel). Each curve is 
an average of a hundred showers. 
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FIG. 9. Angular distribution of the electric field radiated by 1 (short dash), 10 (long dash) and 100 TeV (solid) vertical shower 
crossing from ice to air. The electric field is evaluated at 200 MHz (left panel) and 1 GHz (right panel). Each curve is an 
average of a hundred showers. 


7. Different media 


So far, we have studied the properties of TR from show¬ 
ers crossing an ice-air interface. Salt, soil and regolith are 
also potentially relevant for an experiment, and we per¬ 
formed simulations of 100 TeV showers starting in these 
media. For soil we took parameters of the soil from the 
Argentinian village El Sosneado, next to the Pierre Auger 
Observatory location. The angular distribution of the 
electric field radiated in salt/regolith/soil simulations is 
compared with the standard ice simulation in Fig. |10| 
At 200 MHz, the electric field is comparable for all four 
media. At 1 GHz, a lower emission is obtained in ice. 
This can be explained by the larger Moliere radius of 
this medium (~ roughly double of the other tm), which 
results in a less compact shower. Since the condition of 
coherent emission from the shower particles is then more 


easily broken, a weaker electric field is expected. Soil 
and lunar regolith have rather similar behaviors at both 
frequencies. 


To illustrate an additional feature which can occur 
for certain configurations, we performed a simulation, 
also shown in Fig. in which the refractive index of 
a medium, otherwise similar to ice, was artificially set 
to value of 1.28 so that the Cherenkov angle is smaller 
than the critical angle for total internal reflection in this 
medium (see Sect HI A1) and the coherent Askaryan ra¬ 
diation produced by the shower in the artificial ice is 
now transmitted into air (see Sect HIAl). Notice that 
the refracted Cherenkov radiation is significantly larger 
than the TR, and presents a characteristic angular dis¬ 
tribution, with a beam clearly visible at 1 GHz where 
coherence effects are stronger. A signal of similar origin 
will be encountered in the study of inclined showers in 
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B. Inclined showers 


While some specific applications of transition radiation 
from close-to-vertical showers may be found (for exam¬ 
ple, showers emerging from a mountain range), a more 
likely use of TR could be the detection of UHE neutri¬ 
nos through inclined showers emerging from the Earth 
surface. In fact, the Earth becomes increasingly opaque 
for neutrino energies above 100 TeV [33], and most of the 
observable showers would originate from Earth-skimming 
neutrinos interacting below the Earth surface. 

To gain insight on the properties of TR from inclined 
showers, we performed simulations of 100 TeV showers 
with large zenith angle crossing the boundary from ice 
to air. The shower’s starting point in ice was chosen 
so that the shower maximum development would occur 
at the intersection of the shower axis and the boundary 
plane. Results of the simulations for two showers with 
incident angles of 0sh = 45° and 75°, respectively, are 


first shown in Figs. and 1^ In these graphs, each 
point corresponds to the electric field in a direction (0, 
of observation, the polar angle of the plot being equal to 
(j) and its radial distance from the origin proportional 
to 9. Due to symmetry of the situation, we show only 
half of the sky in each plot. The direction of the shower 
axis (0sh,</'sh = 0°) is marked by a cross. The electric 
field was evaluated at 50 and 1 GHz. In Figs. 13 and 14 
we show the electric field for the same physical situation 
evaluated along lines of constant (p tor (j) = 0°,30°,60° 
and 90°; profiles along p and p + 180° are plotted using 
a single curve with the latter corresponding to negative 
values of 6. Each curve is an average over a hundred 
showers. 

Several distinctive features characterize the radiation. 
At 50 MHz, the electric field is sizable over a large por¬ 
tion of the hemisphere even for the 75° shower, with an 
intensity comparable to that of a vertical shower. Also, 
there is a small angular region around the direction of the 
shower axis where the electric field increases by a factor of 


about three, clearly visible in Fig. IR andpT| This peak 
corresponds to Cherenkov emission in air by high energy 
shower particles moving close to the shower axis (see left 
panel of Fig. [^. As also evident in Figs. [T^and[l4| the 
electric field drops significantly for (0, p) =^0sh, 0^ due 
to vanishing t;_L in Eq. ([^ for V II k. 

At 1 GHz, the emission is overall more beamed to¬ 
wards the direction of the shower axis, and this effect is 
more pronounced the higher the frequency and the inci¬ 
dent angle of the shower. For the 75° shower, |R x E| is 
greater than 10“®'^ V/MHz in ^ 50% of the hemisphere 
at 50 MHz, and only in ^ 20% of the hemisphere at 
1 GHz. Another distinctive feature at high frequency is 
the appearance of an angular band with significantly in¬ 
creased signal. The band is clearly evident at 1 GHz, 
where the electric field is an order of magnitude larger. 


Its origin is explained in terms of coherent Cherenkov 
radiation refracted from ice into air. Indeed, we found 
that most of the radiation in this angular region is com¬ 
ing from below the boundary (“refracted contribution”). 
Further insight is given by a simple model where the 
Cherenkov cone produced by a highly relativistic particle 
moving in ice along the shower axis is refracted into air 
according to Snell’s law. The refracted Cherenkov cone 
matches very well the position and shape of the observed 
band, confirming its association to coherent Cherenkov 
emission by shower particles in ice. 

At 10 MHz, the emitted radiation is basically indistin¬ 
guishable for showers leaving either ice or soil when both 
are crossing the interface at their shower maximum. At 
higher frequencies both electric field profiles show rise of 
the Cherenkov peak and increased beaming towards the 
direction of the shower axis; these features start to ap¬ 
pear at lower frequencies in ice for reasons explained in 
section IHI A 71 


IV. CONCLUSIONS 

We have developed a general method to calculate the 
electric field radiated by the ensemble of particles of a 
high-energy electromagnetic shower developing in two 
different media. The algorithm is based on the stan¬ 
dard ZHS Monte Carlo approach, where a particle track 
is divided into sub-tracks, each contributing to the total 
electric field. In the special case of a particle crossing the 
two media, the particle track is split at the boundary, 
and the electric field from the two sub-tracks is evalu¬ 
ated taking properly into account reflection and refrac¬ 
tion at the boundary plane. Transition radiation natu¬ 
rally arises from this procedure, in addition to the coher¬ 
ent Cherenkov radiation produced by the shower particles 
within each medium. 

Our ZHS-TR algorithm has general applicability and 
could be used, for example, to evaluate the radiation from 
an Extensive Air Shower hitting ground, as seen by an 
observer placed above (or below) the Earth surface. In 
this paper, we focused our studies on the radiation emit¬ 
ted at radio frequencies by showers transitioning from a 
dense medium into air. This configuration is relevant for 
a shower originated by a UHE neutrino interaction below 
the Earth surface, escaping ground into the atmosphere. 
The ZHS Monte Carlo treats only electromagnetic show¬ 
ers, which are appropriate for a charged-current inter¬ 
action of a high-energy electron neutrino. However, we 
expect our results not to change significantly for hadronic 
showers (produced in charged or neutral-current interac¬ 
tions of neutrinos of any flavor) since a large fraction of 
their energy is ultimately dissipated in electromagnetic 

processes Hoi- 

Fundamental properties of the radiation were deter¬ 
mined through Monte Carlo simulations of vertical show¬ 
ers crossing from ice to air. We found that the emis¬ 
sion is fairly isotropic, with an electric field strength of 
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FIG. 10. Angular distribution of the electric field radiated by a 100 TeV vertical shower crossing from ice (blue), salt (red 
dashes), regolith (black short dashes) and soil at El Sosneado (green dot-dashes) to air. Also shown is the electric field for 
a simulation where the ice refractive index was set to the unphysical value of 1.28 (blue long dashes). The electric field is 
evaluated at 200 MHz (left panel) and 1 GHz (right panel). Each curve is an average of a hundred showers. 



FIG. 11. Angular distribution of the electric field radiated by a 
100 TeV inclined shower crossing from ice to air with incident 
zenith angle of 9ah = 45° and <j)sh = 0°. Each point corre¬ 
sponds to the electric field in a direction (6, cj)) of observation, 
the polar angle is equal to (j) and the radial distance from the 
origin proportional to 8. The shower leaves medium 1 (ice) 
into medium 2 (air) while crossing at shower maximum. Due 
to symmetry of the situation we show ony half of the sky in 
each plot. The direction of the shower axis is represented by a 
cross. The electric field is evaluated at a frequency of 50 MHz 
(top) and 1 GHz (bottom). 

^ 5 X 10“^^ V/MHz/m for a 100 TeV shower observed at 
10 km distance. The electric field scales approximately 
linearly with the number of particles crossing the bound¬ 
ary and remains coherent up to ^ 1 GHz. At EeV en¬ 
ergies the growth of the signal magnitude with shower 
energy is slower than linear, due to the LPM effect. Sim¬ 
ilar results were obtained when using salt, soil or regolith 
as a dense medium. We also studied the characteristics of 



FIG. 12. Same as Fig. for a shower with incident zenith 
angle of 8sh = 75° and = 0°. 


radio emission from inclined showers. The radiation pat¬ 
tern is quite broad, but becomes more beamed towards 
the direction of the shower axis for larger incident angles 
or higher emission frequencies. The electric field strength 
is similar to that of a vertical shower, but increases sig¬ 
nificantly at frequencies of ~ 1 GHz in a limited angu¬ 
lar band, corresponding to a coherent Cherenkov cone 
refracted from ice into air. Electric field for showers 
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FIG. 13. Angular distribution of the electric field radiated 
by a 100 TeV inclined shower crossing from ice to air with 
incident zenith angle of 9ah = 45° and (j>ah ~ 0°. Each curve 
corresponds to the shower profile with fixed observer position 
0 while varying 0 (see caption of Fig. 111. Plotted for (j) = 
0°, 30°, 60° and 90°; profiles along <(> and (j) + 180° are plotted 
using a single curve with the latter corresponding to negative 
values of 9. The shower leaves medium 1 (ice) into medium 2 
(air) while crossing at shower maximum. The electric field is 
evaluated at a frequency of 50 MHz (top) and 1 GHz (bottom). 
Each curve is an average of a hundred showers. 


FIG. 14. Same as Fig. for a shower with incident zenith 
angle of 0sh = 75° and 4>Bh = 0°. 


leaving soil show similar features, only their onset ap¬ 
pears at higher frequencies due to better coherence in 
this medium. 

Amongst these results, the angular distribution and 
frequency spectrum of the radiation are particularly in¬ 
teresting. Given the wide solid angle of the emission, a 
large aperture experiment becomes feasible. Also, sub¬ 
stantial radiation in the GHz range facilitates detection, 
thanks to the low radio background and advantages in 
detector design at these frequencies. These observations 
encourage further studies to evaluate the potential of a 
large-aperture UHE neutrino experiment based on detec¬ 
tion of transition radiation. 
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